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Abstract 

Single-molecule localization microscopy of biological samples requires a precise knowledge of the employed fluorescent 
labels. Photoactivation, photoblinking and photobleaching of phototransformable fluorescent proteins influence the data 
acquisition and data processing strategies to be used in (Fluorescence) Photoactivation Localization Microscopy ({F)-PALI\/1), 
notably for reliable molecular counting. As these parameters might depend on the local environment, they should be 
measured in cellulo in biologically relevant experimental conditions. Here, we measured phototransformation quantum 
yields for Dendra2 fused to actin in fixed mammalian cells in typical (Fj-PALIVl experiments. To this aim, we developed a data 
processing strategy based on the clustering optimization procedure proposed by Lee et al (PNAS 109, 17436-17441, 2012). 
Using simulations, we estimated the range of experimental parameters (molecular density, molecular orientation, 
background level, laser power, frametime) adequate for an accurate determination of the phototransformation yields. Under 
illumination at 561 nm in PBS buffer at pH 7.4, the photobleaching yield of Dendra2 fused to actin was measured to be 
(2.5±0.4)xl0~^, whereas the blinking-off yield and thermally-activated blinking-on rate were measured to be 
(2.3±0.2)x10~^ and 11.7±0.5 s~\ respectively. These phototransformation yields differed from those measured in poly- 
vinyl alcohol (PVA) and were strongly affected by addition of the antifading agent 1,4-diazabicyclo[2.2.2]octane (DABCO). In 
the presence of DABCO, the photobleaching yield was reduced 2-fold, the blinking-off yield was decreased more than 3- 
fold, and the blinking-on rate was increased 2-fold. Therefore, DABCO largely improved Dendra2 photostability in fixed 
mammalian cells. These findings are consistent with redox-based bleaching and blinking mechanisms under (F)-PALM 
experimental conditions. Finally, the green-to-red photoconversion quantum yield of Dendra2 was estimated to be 
(1.4±0.6)x10~^ in cellulo under 405 nm illumination. 
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introduction 

Single-molecule localization fluorescence microscopy [1,2,3,4] 
has become a central tool to investigate cells with unprecedented 
spatial resolution. In (Fluorescence) Photoactivation Localization 
Microscopy (F)-PALM (hereafter referred to as PALM), genetically 
encoded phototransformable fluorescent proteins (PTFPs) are used 
to sequentially localize individual emitters with sub-dif&action 
accuracy, so as to reconstruct a super-resolved image of pointiUistic 
character. During data-acquisition, the density of emitting 
fluorescent molecules depends on a number of photophysical 
processes, including photoactivation, reversible blinking and 
irreversible bleaching. A thorough characterization of these 
processes is particularly crucial in quantitative studies aiming at 
counting molecules, as recently shown for proteins forming 
multimeric cell machineries [5,6], found in organelles of sub- 



dilTraction size [7] or arranged in clusters smaller than the 
dilTraction limit [8,9]. (For a review see [10]). 

Li this work, we establish the feasibility of measuring 
phototransformation quantum yields of PTFPs directly in cells 
using standard PALM microscopy. A few previous studies have 
measured photobleaching yields of non-phototransformable fluo- 
rescent proteins at the single-molecule level, typically in vitro to 
ensure a sufficientiy sparse molecular density [11,12,13,14,15]. A 
wide range of values (10 -10 ") were obtained, sometimes for 
the same protein using different protocols [12,13]. It was proposed 
that the experimental conditions and notably the employed single- 
molecule immobilization technique could alter the photophysical 
properties of the FP under study. In the case of green-to-red 
photoconvertible fluorescent proteins (PCFPs), an important 
PTFP's subclass of common use in PALM microscopy, an estimate 
of the photobleaching quantum yield of mEosFP was obtained 
(<I) = 3x10 ■') using BSA .surface-immobilization [16], and it was 
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suggested in the case of mEos2 and pcDronpa2 [1 7] that bulk 
measurements of photobleaching yields may not be valid in single- 
molecule microscopy conditions. The photoconversion quantum- 
yield of PCFPs was measured in the case of KikGR [18] and its 
monomeric variant mKikGR [19] in solution or in gels, giving 
values of the order of ~5xlO ''. Other careful studies on e.g. 
Kaede [20] or EosFP [16] did not provide explicit values for 
photoconversion quantum yields. Recendy, estimates of photo- 
conversion efficiencies, another important parameter accounting 
for the total fraction of fluorescent molecules able to undergo a 
transition to the active (red) state, were obtained in cellulo for 
various PCFPs, but again no quantum yield was provided [21]. 

Since the seminal work by Dickson et al [22], the blinking 
behavior of a number of FPs was investigated by fluorescence 
correlation spectroscopy [12,23,24,25], concentrating on fast 
(~kHz) rather than slow blinking regimes (~ 1-10 Hz). The 
occurrence of long-lasting dark states has been noticed [26,27,28], 
but no quantum yields for the corresponding on-qff transitions have 
been reported. Long-lasting dark states of PTFPs deserve special 
attention because they are detrimental to PALM microscopy and 
single-particle tracking experiments. In PALM, clustering and 
molecular overcounting artifacts induced by slow blinking of 
mEos2 or PA-GFP have been recognized in a number of seminal 
studies [27,29,30]. Indeed, transient dark states with lifetimes 
exceeding the typical frametime in PALM data acquisition (10- 
50 ms) cause single emitters to appear several times during a given 
data acquisition. To correct for these blinking-induced artifacts, 
several strategies have been proposed based on grouping 
fluorescence bursts that appear close in space and time. Annibale 
et al [3 1] used an approach in which the fluorescence traces from 
spatially close spots are subdivided whenever an off-time exceeding 
a given threshold tj is reached, resulting in a number of sub-traces 
N(tj). The true number of molecules is then obtained by fitting the 
J^(tj) curve with a semiempirical model accounting for blinking. A 
variant of this technique based on grouping spots with Kalman 
filtering was developed by Lando et al [32]. To accurately group 
spots in both the temporal and spatial dimensions, Coltharp et al 
[33] used another empirical strategy based on calculating a 
Jaccard index as a function of two thresholding variables, txi,„, and 
driims- The optimum t-fhres ^nd d-fhns arc tlu'n applied to the entire 
data set, independent of the local molecular density. Recently, an 
elegant method was introduced by Lee et al [8] that proposes to 
subdivide molecular fluorescence traces based on a cutoff value 
which is dependent on the FP photobleaching and photoblinking 
rates as well as the local density of molecules. The optimum cutoff 
value z,. „i,i giving rise to a minimal counting error is obtained by an 
iterative procedure, the actual number of molecules in each 
evaluated spot cluster being the unknown. In this approach, 
bleaching and blinking rates of the PTFPs are obtained from in 
vitro experiments in which the molecular density can be kept sparse 
enough to prevent spatial overlap between molecules. None of the 
above techniques, however, was designed to extract the actual 
photoblinking quantum yields of the employed PTFPs in cellulo. 

Overall, although the importance of evaluating PTFP's 
photophysics in biologically relevant experimental conditions is 
increasingly recognized [7,32,33], a method to quantitatively 
evaluate their phototransformation yields in densely labeled cells is 
stiU lacking. Such in cellulo evaluation would facilitate accurate 
molecular counting and have a fundamental interest for the 
mechanistic investigation of FP photophysics and the future design 
of improved variants. 

On the structural standpoint, the photoconversion mechanism 
of green-to-red PCFPs such as Kaede [34], Dendra2 [35,36], 
mEos2 [37], mKikGR [19], mClavGR2 [38], mMaple [39], 



pcDronpa2 [17] and their variants has been well characterized, 
involving a light-induced breakage of the protein backbone next to 
the chromophore which results in irreversible elongation of its 

electron conjugation through a fi-elimination reaction 
[16,34,40,41]. In contrast, the precise blinking and bleaching 
mechanisms of these FPs remain incompletely understood, despite 
recent structural studies of IrisFP, a variant of EosFP [42,43,44]. 
Overall, investigations of blinking in FPs has been lagging behind 
those of organic dyes typically used in (direct) STochastic Optical 
Reconstruction Microscopy ((d)STORM) [45,46,47]. A large body 
of work has demonstrated that organic dye photophysics is 
strongly dependent on environmental parameters such as pH, 
presence of molecular oxygen or redox partners, attachment to 
protective agents [45,47,48] and even the nature of the molecule 
labelled by the dye [49]. Although the 4-(/»-hydroxybenzylidene)-5- 
imidazolinone (/)-HBI) chromophore in fluorescent proteins is to 
some extent protected from the external medium by the fi-barrel 
matrix, recent evidence has been provided that such environmen- 
tal dependence also applies to FPs. Redox agents [50,51,52], pH 
[23,25], temperature [53], physicochemical parameters such as 
viscosity [54] and even external electric or acoustic fields [55] 
indeed have been shown to greatly modulate the photoconversion, 
photoblinking, photoswitching or photobleaching behavior of a 
number of FPs. Likewise, fixation conditions or type of fusion 
construct employed in biological cells could play a stringent role 
[56]. 

Performing in cellulo single-molecule measurements of PTFP's 
photophysics is not straightforward. Such measurements are 
typically hindered by unfavorable signal-to-noise ratio of the 
detected fluorescence, and by the potentially high density of 
labeled molecules. In PALM, due to residual activation by the 
readout laser beam, a sufficiently sparse distribution of activated 
molecules in each frame is not guaranteed. To extract photo- 
transformation yields, the photoexcitation rates of individual 
molecules need to be assessed, which requires the precise 
knowledge of not only their extinction coefficients but also of the 
local laser power-density and polarization (which in turn 
complicates the use of a TIRF geometry). Fluorescent molecules 
might adopt a variety of orientations, and the strong dependences 
on dipole orientation of both photon absorption and fluorescence 
response need to be taken into account. Finally, the inhomogeneity 
of the chemical nanoenvironment in the cell may result in multiple 
photophysical behaviors coexisting within the sample. 

Here, we propose a method to tackle these difficulties. To 
estimate the phototransformation quantum yields of green-to-red 
PCFPs in cellulo, we devised a strategy based on an extension of the 
%c optimization approach of Lee et al [8] . The performance of our 
approach was evaluated by performing extensive simulations, so as 
to establish a range of experimental conditions for which the 
method is expected to work reliably. The phototransformation 
yields of Dendra2 fused to actin were then measured in fixed 
mammalian cells, and compared to those of Dendra2 embedded in 
a poly-vinyl alcohol (PVA) matrix, in the presence or absence of 
the antifading agent l,4-diazabicyclo[2.2.2]octane (DABCO). 

Methodology 

Generation of simulated data 

To generate simulated data, comprehensive PALM simulation 
software was written in Matiab (Figure SI, Scripts SI), with 
additional features as compared to similar software developed to 
date (see e.g. [31,33,57]). A chosen number of dye molecules are 
randomly placed on a defined image pattern (typically a 
segmented nanoscopy image), so as to control dye molecular 
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Figure 1. Photophysical model of Dendra2 phototransformations. The inactive state represents the neutral state of the green chromophore 
whereas the active state represents the anionic state of the photoconverted red chromophore. Fast exchange is assumed between protonated and 
neutral states of the chromophore. Possible bleaching pathways through the reversible dark state are neglected, as well as blinking and 
photobleaching pathways involving inactive molecules (in gray). 
doi:1 0.1 371/journal.pone.0098362.g001 



density. The photophysical properties of the fluorescent markers 
(mostly phototransformable fluorescent proteins) are defined in a 
database containing parameters such as excitation and fluores- 
cence emission spectra, fluorescence quantum yield, phototrans- 
formation action spectra as well as photoconversion, on-off 
photoswitching, photobhnking and photobleaching quantum 
yields. Properties of the PTFP's chromophores in their anionic 
and neutral states are distinguished and simulations at various pHs 
can be performed. The orientation of the chromophore dipoles 
can be selected as "tumbling" (emitters are freely rotating, 
representative of e.g. a live cell), "frxed" (emitters assume a 
random orientation, but constant over time, representative of e.g. 
a fixed cell), or "uniform" (emitters adopt the same dipole 
orientation, constant over time, representative of e.g. oriented 
molecules in a fixed or live cell). 

The experimental setup is simulated by a number of parameters 
describing laser illumination, microscope optical properties and 
detector characteristics. Four laser beams with Gaussian spatial 
profiles within the focal plane can be defined (i.e. for activation, 
switching, back-switching and readout), with chosen wavelengths, 
powers, polarization (circular or linear) and FWHM (full width at 
half maximum). Complex laser sequence patterns can be 
generated during data collection, such as for example pre- 
bleaching steps with the readout laser, and ramping-up or a 
Fermi profile [8] for the activation laser. The microscope is 
described by an objective of given numerical aperture JVA, which is 
assumed to give rise to a Gaussian point spread function (PSF) with 
standard full width at half maximum {FWHM= 1.22)JMA where X 
is the emission wavelength). The photon collection efficiency takes 
into account the orientation of the emitting dipoles based on 
Fourkas et al [58]. Transmission efficiencies of inserted emission 
filters, as well as overall transmission efficiency of the optical setup 
are also taken into account. The EMCCD detector is described by 
the effective pixel size, dark current noise, readout noise and gain 
(counts/detected photon). 

PALM data sets with a defined frametime and number of 
frames can then be generated. For each single-molecule emission, 
the spatial (throughout the microscope PSF) and spectral 
(throughout the emission spectrum) distributions of emitted 
photons are calculated based on poissonian statistics. A defined 
level of poissonian autofluorescence background, optionally 
proportional to the local illumination power density, can be input. 
The program outputs the stack of acquisition frames, diffraction 
limited images, laser illumination profiles, and for each activated 
molecule, its x andj coordinates, fluorescence emission trace, the 
number of detected photons and an estimate of the signal-to-noise 
ratio. 



Additional features, not of direct use for the current work, 
include possibilities to simulate dual-channel detection, two-color 
experiments (with two different fluorescent markers), sample drift, 
single-molecule Forster resonance energy transfer (FRET) between 
two fluorescent markers, and local activation/bleaching similar to 
what can be achieved with a Fluorescence Recovery After 
Photobleaching (FRAP) module. 

Biophysical model 

To extract the phototransformation yields of green-to-red 
photoconvertible fluorescent proteins, the photophysical model 
of Figure 1 is assumed. The inactive state of the marker 
corresponds to the green protonated chromophore whereas the 
active state corresponds to the red anionic chromophore. The 
green anionic chromophore (non-photoconvertible) and the red 
neutral chromophore (non-fluorescent) are thought not to play a 
role and are therefore not considered. However, rapid exchange 
between anionic and neutral chromophores relative to the 
acquisition frametime is assumed so that the rates of excitation 
in either the inactive or active states are evaluated based on 
extinction coefficients experimentally measured at a given pH. 

Our photophysical model assumes that pathways leading to 
reversible dark states and to nonreversible bleached states are 
independent, similar to the approach taken by other groups [8,33]. 
In doing so, we thus neglect the hypothesis that photobleaching 
may follow from photochemical transformation of a long-lived 
dark state. Although such an hypothesis has been proposed in 
several photophysical schemes [59,60], we adopt this simplification 
based on our recent mechanistic investigations of IrisFP for which 
the major photobleaching pathways were suggested to differ from 
that leading to long-lived dark states [44]. 

In the model of Figure 1 , the recovery from the dark to the 
active state is characterized by a single rate k^ji and does not 
incorporate a quantum yield. This stems from the fact that in the 
absence of 405 nm light (readout activation mode) or in the 
presence of weak 405 nm light (typically used in PALM data 
acquisition), the process is essentially thermally induced (ref [8] 
and our own data (not shown)). 

Finally, the model also assumes that blinking and bleaching 
following excitation of the inactive chromophore are negligible 
(Figure 1, left gray part). Such an assumption is justified because 
typical photoconversion yields measured to date are at least one 
order of magnitude higher than blinking and bleaching yields. 
Thus, the pool of inactive molecules competent for photoconver- 
sion is expected to be only marginally reduced at a given time. In 
addition, the various states of the inactive chromophore cannot be 
easily separated, as none of them fluoresces upon excitation by the 
readout laser. 
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Figure 2. Sketch of the processing pipeline. 

doi:10.1371/journal.pone.0098362.g002 



Experimental setup 

A home-built PALM set-up was used based on an Olympus 
1X81 microscope (Olympus, Japan) equipped with diode-pumped 
solid-state 405 nm (CrystaLaser LC, USA), 488 nm (SpectraPhy- 
sics, USA) and 561 nm (Cobolt Jena, Sweden) lasers. Widefield 
illumination was achieved by focusing the laser beams to the back 
focal plane of a 60 x 1 .45 MA oil immersion apochromatic objective 
lens (Olympus, Japan). The intensity and time-sequence of laser 
illumination at the sample was tuned by an acousto-optical tunable 
filter (AOTF, Qrianta Tech, USA). Near circular polarization of 
the laser beams was ensured by inserting a polychromatic quarter- 
wave plate downstream the AOTF. To prevent sample drift 
during data acquisition, the samples were placed on an IX2-NPS 
"nosepiece" stage (Olympus) frxed directly on the objective. 
Fluorescence images were acquired with an Evolve back-illumi- 
nated EMCCD camera (Photometries, USA) controlled by the 
Metamorph software (Molecular Devices, USA). 

Sample preparation 

In vitro samples. Purified Dendra2 was diluted to nanomo- 
lar concentration in a pre-photobleached 1 % solution of PVA in 
PBS at pH 7.4. Where appropriate, l,4-diazabicyclo[2.2.2]octane 
(DABCO) (Sigma, USA) was added to a final concentration of 
10 mM. The solution was spread as a thin film on a coverslip and 
let to harden at ambient temperature before PALM data 
acquisition. Control experiments with PVA devoid of Dendra2 
showed no blinking single molecule (Figure S2). 



In cellulo samples. Mammalian cells were grown on 
conventional coverslips or LabTek II 8-well chambered coverslips 
^1.5 (Nalgen-Nunc, USA). Before seeding the cells, LabTek 
chambers or coverslips were washed 3 times with 80% acetone in 
water, 100% methanol and milliQ_ water, immersed in 1 M 
solution of KOH overnight, and thoroughly washed with miUiQ, 
water. For expression of Dendra2-actin chimera in mammalian 
cells (HeLa Kyoto and Vero), we used the plasmid Dendra2-actin/ 
pMCl (described earlier [61]), which encodes for a humanized 
version of the Dendra2 gene from the pDendra2-C vector 
(Evrogen, Moscow, Russia). Cells were grown in complete 
Dulbecco's modified Eagle's medium (DMEM, Life technologies, 
USA) (4.5 g/L glucose, no phenol red) supplemented with 10% 
heat-inactivated fetal bovine serum (Life technologies), 2 mM 
glutamine, 100 Hg/ml penicillin, 100 U/ml streptomycin and 
1 mM sodium pyruvate (all from Life technologies), in 5% CO2 
atmosphere at 37°C. No vitamin was added to the cell culture 
medium, which however already contains a complete set of 
vitamins, including riboflavin and pyridoxal [52,62]. The cells 
were transfected with Dendra2-actin/pMCl plasmids 
(0.5 |a.g/well) using Lipofectamine2000 (Invitrogen, USA) or 
FuGeneHD (Roche, Switzerland) as transfection reagents. At 24 
to 48 hours post-transfection, the ceUs were washed with PBS, 
fixed with 4% paraformaldehyde (Sigma, USA) at 20°C for 30- 
60 min, washed with PBS 3 times, permeabUised with a 0. 1 % 
solution of Triton x 100 in PBS, and washed with PBS for another 
3 times. Where appropriate, DABCO (Sigma, USA) was added to 
a final concentration of 10 mM. CeUs were kept under pre- 
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photobleached PBS at pH 7.4 during PALM data acquisition. The 
molecular oxygen level was in equilibrium with ambient atmo- 
sphere. The coversUps were mounted on microscopy slides using 
Parafilm gaskets and sealed with nail polish. 

Data acquisition 

Reference laser beam profiles were recorded before each 
experiment using a coverslip uniformly marked with a fluorescent 
pencil. PALM datasets were acquired at room temperature (20- 
22°C) using two illumination schemes: 1) continuous illumination 
with 561 nm laser, which served both to photoconvert Dendra2 
fluorophores and to readout the fluorescence of the activated red 
form; 2) continuous illumination with 561 nm and 405 nm laser. 
The power density at the focal plane was measured in each 
experiment and was typically ~5-7 kW/cm^ at beam center for 
the 561 nm laser and ~0.1-1 W/cm^ for the 405 nm laser. Data 
sets were acquired at 30 ms/frame. 

Processing pipeline 

To determine the on-off blinking quantum yield, the off-on 
blinking rate, the photobleaching quantum yield as well as the 
photoconversion "brightness" (see below) of the fluorescent 
marker, the processing pipeline sketched in Figure 2 was 
developed in Madab (Scripts SI). 

Single-molecule localization is first performed based on 
standard 2-D Gaussian fitting of the diffraction-limited fluores- 
cence spots, after rolling ball background subtraction. In addition 
to the precise localization of the fluorescent molecules, attention 
here is focused on extracting their accurate fluorescence traces. 
Indeed, removing a poorly localized spot from the dataset may 
result in introducing a glitch in a fluorescence trace, which would 
later be interpreted as a blinking event. Such problems may occur 
in the case of low signal-to-noise ratio or when multiple spots 
overlap at high molecular density. We modified our PALM 
localization software to carefully select spots judged appropriate 
for trace reconstruction (e.g. by setting criteria to allow center-of- 
mass based localization instead of rejection in the case of poor 
Gaussian fitting). 

The next step consists of grouping localized spots found in 
subsequent frames, so as to constitute initial clustc-rs of spots from 
which the most probably correct number of individual molecules 
and their associated fluorescence time traces can be extracted 
based on splitting. To achieve initial clustering, localized spots 
are first connected in space and time according to coarse threshold 
criteria d(.si,.S2)<5a, where d(si,S2) is the distance between spots Si 
and $2 and a is the localization uncertainty, and t{si,S2)<5x where 
t(si,S2) is the time separation between spots .sy and S2 and T is the 
maximum expected molecular lifetime. (We noticed that a more 
stringent distance criteria (e.g. Scr) repeatedly leads to misclustering 
events.) However, such connection patterns can quickly lead to 
large groups of spots that spread over unrealisticaUy large 
distances, especially in densely labeled and continuous biological 
structures. To subdivide such groups of spots into meaningful 
initial clusters, we developed an algorithm based on kmeans 
clustering analysis. Figure S3 shows an example of how densely 
packed localized spots are distributed into initial clusters. 

Once initial clusters are determined, their fluorescence time 
traces can be built. We then use an approach similar to that 
introduced by Lee et al [8] to optimally estimate the number of 
molecules that constitute each cluster and extract their respective 
time traces. This approach consists of splitting cluster traces at 
every ofl-time T,^that is larger than a cutoff value T^, so that the 
molecular counting error is minimized. The optimum value has 
been shown to depend both on the blinking and bleaching rates of 



the fluorescent label and on the number jV of labels present in the 
cluster [8]. As jVis unknown, an iterative procedure was proposed, 
based on the computation of an optimal "tc-curve", T^ ,,^,, ^K-N), 
which was shown to rapidly converge to the correct number of 
molecules. Importantiy, Lee et al construct their T,-curve based on 
photophysical parameters of the fluorescent markers that are 
independendy measured in vitro. In marked contrast, our goal is 
precisely to assess these parameters in celluh. Thus, the following 
global iterative procedure was developed (Figure 2, orange dashed 
line). First, an initial fixed Xc value (e.g. 0.2 s) is used to split all 
cluster traces in the dataset so as to obtain a starting set of 
individual molecular fluorescence traces. These individual traces 
are then processed to extract bleaching and blinking histograms: 
■Nbkack=h(T,,i^^, N„n-h(TJ, Noff=h(T„J, where Tj^^^ r„„ and 
represent the time before irreversible bleaching, blinking off, 
and blinking on, respectively, and .Ni,k„ch, J^<m and j\v^ represent the 
number of events within given ranges of T^^^^/,, T„„ and 7)^. Based 
on our photophysical model (Figure 1), and on the proper 
evaluation of the spatially- and dipole-orientation-dependent 
molecular excitation rates (see below), these histograms can be 
fitted to extract a first evaluation of bleaching and on-o^ blinking 
quantum yields, as well as a rate for thermally activated off-on 
blinking. The self-consistency between these quantum yields can 
be verified by computing a blinking probability histogram, 
P = h(Mui„^), where Muint represents the number of blinking events 
per molecule, and which is expected to follow a geometric 

distribution P{Nhnnk = n)=r]"{\-ri) with )?=- [8]. 

'VdA + I'm 

From the knowledge of <1>m, ^,la, and 4^4 it is then possible to 
generate a first T^-curve, which is used to split again the cluster 
traces. The whole procedure is repeated until convergence, which 
typically requires less than 5 cycles. At the end of the process, final 
quantum yields and associated standard deviations are obtained, 
together with a quantitative evaluation of the total number of 
activated molecules. The number of newly activated molecules in 
each frame also becomes available, which allows for an estimation 
of the photoconversion brightness (sec below) and the total 
number of (correctiy matured) labeling molecules in the sample. 

Evaluation of excitation rates 

To accurately fit bleaching and blinking histograms, a proper 

evaluation of the spatially- and dipole-orientation-dependent 

molecular excitation rate is necessary. The excitation rate [s '] 

of a rapidly tumbling molecule located at position (x^) is given by: 

, (IO-'^)Lm(IO) _, _, 

k = eP/.- where s [M cm 1 is the extinction 

NaHc 

coefficient of the considered species (/ or A molc-cules) at the 
excitation wavelength X [nm], P [Wcm"^] is the laser power- 
density at (x^), Na is the Avogadro number, h is the Planck 
constant and c is the speed of light. In this work, chromophore 
extinction coefficients were derived from published reference 
values and fluorescence excitation or absorbance spectra at given 
pHs [61]. 

For a fixed oriented molecule with spherical coordinates 6 and 

<p, the excitation rate of the absorbance dipole is angle dependent. 

For a circularly polarized laser beam whose electromagnetic field 

is assumed to be parallel to the objective focal plane, the excitation 
3 

rate is given by: kg = -kcos^O, where 9 is measured from the focal 

plane. For a linearly polarized laser beam whose electromagnetic 
field is in the objective focal plane and aligned with the x axis, the 
excitation rate becomes: ko,p = 'ikcos^9cos^(p (see Appendix SI). 

To handle the spatial heterogeneity of the laser beam under 
widefield illumination conditions, lasers operating in the TEMqo 
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mode are assumed. In such a case, the intensity distribution is 
Gaussian, and experimentally recorded laser beam profiles can be 
fitted for their center x„ and their 7// radii r„ r^,. The power 
density at every single molecule positioned at xo,yo (taken as the 
average of the x, y coordinates of the localized spots constituting 
the molecule) is then obtained as: 



2Pt 

P(a-o,1'o)= exp 



where Pf is the total power measured at sample position. We note 
that the assumption of Gaussian beams is typically adequate in 
widefield illumination conditions, but not in Total Internal 
Reflection Fluorescence (TIRF) mode. In addition, the exponen- 
tial decay of the electric field amplitude along the axial direction in 
TIRF mode and its complex polarization pattern would prevent a 
reliable evaluation of local excitation rates [63]. As a consequence, 
proper estimation of phototransformation yields is best achieved 
under widefield conditions, although at the expense of a lower 
signal-to-noise ratio. 

In practice, to obtain rehable fits of the Nhbadi and N„„ 
histograms, molecules are sorted in a limited number of annular 
concentric subregions Ri, R2, ... R„ in which P(xy) is considered 
homogeneous with average power density Pj = (^P(x,y)'}j^.. 
Global fits of the photobleaching and on-oflF blinking quantum 
yields are then performed according to the following formula; 



excitation wavelength (typically 561 nm), £5^; (Figure S3). This 
stems from the fact that, in this wavelength region, the absorbance 
spectrum of protonated green molecules (which exhibits a broad 
band centered at ~400 nm) cannot be reliably measured, even at 
very low pH, because not only is it very weak but also it is always 
obscured by residual overlapping absorbance from anionic 
chromophores. One solution to get around this problem is to 
report the photoconversion brightness, instead of the photocon- 
version quantum yield, defined as (Pa = E.^Bi'l'a- Indeed, the quantity 
that can be accurately fitted is k/^^, and thus the photoconversion 
brightness is given by ^a = 5.21 xlO*A/Da/P/,. Overall, similarly to 
the fitting procedure for retrieving photobleaching and on-oflF 
blinking quantum yields, the photoconversion brightness and total 
number of labeled molecules are obtained according to the 
following formula: 

(') = 'Vj X [1 - exp( - 1 .92 X 10~' X //>'■ X f)„/)] (lumhling) 
!^'a(i)=K>^ I ' coso|l-exp^-^cos^0xl.92xl0-^x/.i>'x(!)„/jja'0 
(fixed, circular) 

fli fi/2 

N'Jl)=N'„x\ cos0[l-exp(-3cos^0cos^i?)xl.92xl0-*xA/"x5i„/)] 
Jo Jo 

(fixed Jinear) 



Ni(t) = N^ki<l>e\p{~kiOr) (tumbling) 

^,(0 = ^A'o'w3'x J ' cos^ O£xp(^-^kicos^0x<i>?jd0 

{fixed, circular) 

p2n pn/2 

Ni{t) = 3Nl^ki(S>x cos-' Ocos^ipe\p{-3kicos^0cos^(px<l>t)d()d(p 

Jo Jo 

[fixed, linear) 

Obtaining the thermally activated rate for return from the dark 

state, krA, is simpler because there is no dependence on 
chromophore position or dipole orientation: assuming a homoge- 
neous population of molecules: N{t) = NokrASxp{ — krAt). 

Determination of photoconversion brightness under 
readout laser illumination 

As already mentioned, once optimal splitting of the cluster's 
traces is achieved, the number of newly activated molecules per 
frame along data collection can be extracted. In the presence of 
readout excitation only, and assuming constant illumination 
during the acquisition time, the green-to-red photoconversion 
yield (I>„ as well as the total number of labeled molecules present in 
the investigated sample area can in principle be obtained. In the 
case of rapidly tumbling molecules, the accumulated number of 
activated molecules along time MaA) can be fitted according to: 

NA{t) = No[\-exp{-k,Q>j)] 

where No is the total number of labeled molecules and kj is the 
excitation rate of inactive molecules by the readout beam. 
However, a difficulty is that it is practically impossible to 
accurately estimate the extinction coefficient of the inactive 
molecules (green protonated chromophores) at the readout 



Determination of photoconversion quantum yield under 
405 nm activation laser 

In the presence of weak 405 nm activation, the same procedure 
as described above can be used as long as the density of activated 




10'' 

10' 10^ 10^ 10"* 

Input blinking-off yield 



Figure 3. Accuracy of bllnking-off quantum yield retrieval. 

Extracted blinking-off quantum yields are plotted against a range of 
input yields used in a series of simulated PALM data sets. Those 
simulated data sets were generated using parameters displayed in 
Table SI (except for the varying blinking off input yield). Error bars 
represent fitting errors. Yield retrieval for highly blinking molecules 
underestimates the true values, due to the limited time-resolution in 
PALM experiments. For weakly blinking molecules, retrieval accuracy is 
superior, but precision is worsened (larger error bars) because of 
statistical uncertainty due to low numbers of blinking molecules. The 
inset shows the same curve in an inverted manner. A polynomial fit to 
the data points provides a potential calibration curve allowing, in real 
experiments, to estimate the true yield from the extracted yield. 
doi:10.1371/journal.pone.0098362.g003 



PLOS ONE I www.plosone.org 



6 



June 2014 | Volume 9 | Issue 6 | e98362 



Phototransformation Yields in Localization Microscopy 



— Bleaching yield 
—»— Blinking-off yield 
- BiinKing-on rate 
Photoconversion yield 




5 10 15 20 25 30 
Mean signal to noise ratio per molecule 



1.0 



0.5 



'— 0.0 



-0.5 



D 


1— I 













-1.0 

0 2 4 6 8 10 
■Molecular density (active molecule per frame / pm') 



1 ■ 1 ' 1 ' 1 ' 1 ' 1 


0.5 - 




0.0 - 




-0.5 - 




-1.0,, 
-1.26 - 



10 20 30 40 50 60 
Frame time (ms) 




1000 10000 
Laser power density (W.cm'^) 




0.0 0.2 0.4 0.6 0.8 1.0 
Fraction of photoconverted molecules 



1.0 



0.5 



0.0 



-0.5 



-1.0 



0 3000 6000 9000 

Total number of photoconverted molecules 



Figure 4. Effect of experimental parameters on accuracy of phototransformation yield's retrieval. In each panel, extracted 
phototransformation yields are plotted against a key experimental parameter (signal to noise ratio per molecule (A), frametime (B), laser power (C), 
molecular density per frame (D), fraction of photoconverted molecules (E) and total number of photoconverted molecules (F)). In each case, PALM 
datasets were simulated using parameters of Table SI (except for the varying parameter). Black: bleaching yield, blue: blinking-off yield, gray: 
blinking-on rate, red: photoconversion yield. Yield retrieval accuracy is estimated by the logarithm of the ratio between true and extracted yields and 
error bars represent fitting errors (with 95% confidence). In D, the insets show representative simulated frames at molecular densities of 1 and 8 
active molecules per frame per square micron, respectively. In E, only the extracted photoconversion yield is shown, as the input parameter is not 
relevant to bleaching and blinking. In F, conversely, only extracted bleaching and blinking yields are represented. Experimental data discussed in this 
work have been obtained within ranges of conditions indicated by braces and/or black arrows. Molecular densities per frame cannot be easily 
assessed for experimental data and therefore reported values are only indicative. 
doi:1 0.1 371 /journal.pone.0098362.g004 



molecules in each frame can be managed. An apparent 
photoconversion brightness cpa,app can then be extracted (which 
results from photoconversion by both activation and readout 
lasers), from which the photoconversion quantum yield at 405 nm 

. , . , . 'Pa,ann . Psilhil 

•Pa 405 can be estimated as: Oa,405 « l.Jsx— x 

£405 "405^05 

"'^^^ , where Psei and P^os refer to the power densities at beam 
£405 

center, (567 and Z705 refer to exposure times and 1.38 is the ratio of 
the two wavelengths. 



Results 

Simulated data 

To test the proposed strategy for estimating phototransforma- 
tion yields, we first simulated a representative PALM data set using 
our Matlab-based simulation software. 30000 Dendra2 molecules 
were randomly spread over a microtubule-mimicking pattern 
(corresponding to a molecular density of 545 molecules/ |J.m^) and 
6000 frames with 30 ms exposure were generated using circularly 
polarized 561 nm laser excitation. No 405 nm activation laser was 
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used and, under readout laser activation, the green-to-red 
photoconversion brightness was set to 4.3xl0~* M~'cm~' 
(<5<,= 1.2x10"*; e5ff/ = 3.57 M"'cm~'), so that the average 
density of activated molecules per frame amounted to 0.22/|J.m^. 
The simulated sample was assumed to be representative of a live 
cell, so that each Dendra2 molecuk- was modeled as freely 
tumbling. Input values for the bleaching yield, the blinking-olf 
yield and the thermally activated blinking-on rate were taken as 
2xlO~^, 5xl0~^ and 20 s~' respectively, close to the experi- 
mental values that can be extracted from Lee et al [8] . A complete 
Ust of parameters used to generate this PALM dataset is shown in 
Table SI. The simulated sample is shown in Figure SI. 

The dataset was then submitted to the processing pipeline 
described above. Erroneous phototransformation yields and 
blinking-on rate were deliberately input to the fitting routines as 
starting values (with deviations by an order of magnitude relative 
to true values). Five concentric regions where the laser beam can 
be considered homogeneous were defined as shown in Figure S5. 
Convergence was achieved after 3 cycles and the final retrieved 
values are reported in Table S2 and compared with the "true" 
values used in the simulation. A representative T^-curve is shown in 
Figure S6 and histograms of the correlation coefficients between 
true and retrieved single-molecule traces are presented in Figure 
S7. 

The photobleaching yield was retrieved with high-fidelity, with 
an error less than 1%. Recovery of the blink-on rate was also 
satisfactory (21.9 s 9% error). However, a substantial discrep- 
ancy was obtained between the retrieved and the true blink-olf 
yields. The retrieved yield (3.6x10 **) was smaller than the true 
yield, with an error of 28%. This finding is explained by the 
limited time-resolution of PALM data acquisition (see Discussion 
SI, Figure S8, and see below). Retrieval of the photoconversion 
brightness was fully satisfactory (4% error), allowing at the same 
time to extract an accurate estimation of the number of labeled 
molecules pres(;nt in the sample (8"/!) error), a crucial aspect in 
quantitative super-resolution microscopy [9], although this 
retrieved number does not take into account additional effects 
such as reporter misfolding or presence of a non-photoconvertible 
fraction of the labels [7,21]. 

Effect of dipole orientations and laser polarization 

Two additional mock data sets were generated using the same 
set of parameters except that this time the simulated sample was 
assumed to be representative of a fixed cell: each Dendra2 
molecule was assigned a random dipole orientation, constant ()\ er 
time. In the first dataset a circular polarization of the laser beam 
was assumed whereas a linear polarization was chosen in the 
second one. In both cases, using the appropriate fitting models, the 
quality of yield extraction followed the same trend as for tumbling 
molecules, although the overall accuracy was slightly lower for 
circular polarization and significantly lower for linear polarization 
(Table S3). This deterioration follows from the strong fluctuation 
in the ability of each molecule to absorb photons. The results were 
further degraded if a fitting model corresponding to tumbling 
molecules was used. This can be understood by inspection of 
Figure S9: in the case of molecules with fixed orientations and a 
circularly polarized laser, a few molecules experience a higher 
excitation rate than the average rate experienced by a tumbling 
molecule (i.e. due to an absorbance dipole parallel to the objective 
focal plane) or a lower rate (i.e. due to an absorbance dipole along 
the optical axis). As a consequence, these molecules systematically 
exhibit shorter or longer fluorescence time traces, slighdy 
distorting the shape of the Ni,ieaA and jV„„ histograms. In the case 
of a linearly polarized laser, this distortion is more pronounced 



(Figure S9). Likewise, the shape of the cumulative plot of the 
number of photoconverted molecules along time is distorted when 
molecules adopt fixed orientations, as compared to tumbling 
molecules, and more so when the excitation laser is linearly 
polarized (Figure S9). 

These simulated data highlight the importance of taking into 
account fluorophore dipole orientations in such studies. Never- 
theless, based on the results presented in Table S3, we anticipate 
that, at least when circularly polarized lasers are used, errors made 
by assuming a somewhat incorrect molecular orientation model (as 
could be the case in biological samples) should be relatively small 
in comparison with other sources of errors (see below). 

Retrieval of blinking-off yield 

The noticeable underestimation of the blinking-off quantum 

yield in the above simulations calls for an investigation of the range 

of blinking yields that can be adequately retrieved under typical 

PALM data acquisition parameters. To this aim, we performed a 

number of simulations in which the blinking-off quantum yield 

was systematically varied, assuming a fixed cell with circularly 

polarized laser, and keeping all other parameters identical to those 

described above. The results (Figure 3) show that retrieval of the 

blinking-off yield progressively deviates from (underestimates) the 

true value as the latter is increased. This is explained by the fact 

that, as the input yield is increased, more and more short on-times 

events are missed due to the limited time-resolution of the 

experiment (see Discussion SI and Figure S8). On the other hand, 

-6 ®<W 
for low yields on the order of 1.0x10 (or rather — — <<1), a 

significant statisticEil uncertainty is observed, due to the fact that 
very few molecules experience blinking, thus giving rise to noisy 
jV„„ histograms. Interestingly, Figure 3 suggests that in the range 
~2.0xlO~*'- 1.0x10"* an accurate determination of the correct 
yield could be retrieved if a proper calibration is done. Calibration 
curves of the type shown in Figure 3 may be obtained from 
simulated data in which experimental parameters are reproduced 
faithfully, notably the rate of excitation of the fluorescent 
molecules and the acquisition frametime. Finally, Figure S 1 1 
shows that the correct extraction of the bleaching quantum yield, 
photoconversion brightness and ofF-on blinking rate arc not 
significandy affected by the value of the blinking-ofiF yield, except 
at very high blinking yield > 1 .0 X 1 0 *. 

Evaluation of suitable experimental parameters 

In order to estimate experimental conditions suitable for 
accurate determination of phototransformation yields, a number 
of additional simulations were performed in which the same set of 
yields (representative of Dendra2) were input and the following key 
parameters were systematically varied: background level, frame- 
time, laser power, molecular density and number of frames. The 
results are presented in Figure 4. 

Background leveL High autofluorescence or background 
noise due to out-of-focus light is a serious concern when cellular 
samples are investigated using a widefield geometry. Such noise 
may prevent the detection of many fluorescent spots, which results 
in a significant number of single-molecule traces being chopped 
into shorter traces. This in turn translates into increased apparent 
bleaching and blinking yields. At medium signal-to-noise ratio 
[SMR, defined as the photon output of individual spots divided by 
the standard deviation of the background over the spot area), this 
effect actually compensates for the underestimation resulting from 
time truncation (see above), explaining why errors are not at their 
minimum in the absence of background noise (Figure 4A). 
Importantly, in vitro data reported in this work typically match 
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Figure 5. Representative Kiistograms of Nbieacir N, 
the absence and in the presence of 10 mM DABCO; 
DABCO. 

doi:1 0.1 371 /journal.pone.0098362.g005 



on Noffitom in vitro and in vivo Dendra2 pKiotophysics. A-C: PVA-embedded Dendra2 in 
D-F: Fixed Vero cells expressing Dendra2-actin, in the absence and in the presence of 10 mM 



this medium signal-to-noise ratio [SNR ~15), wliereas in cellulo data 
are sliglitly worse [SNR ~ 10), in a regime where phototransforma- 
tion yields could be slightly overestimated. 

Frametime. reducing the frametime enables enhanced time 
resolution at the expense of reduced signal-to-noise ratio. 
Inspection of Figure 4B reveals that, given the input value we 
have chosen for the yields, a good compromise is maintained in 
the range 10-30 ms frametime. At the shortest frametime, the 
reduced SNR results in a slight overestimation of the photoblinking 
regime, whereas at the longest frametime (60 ms) most fluores- 
cence time traces are contained within a single frame, completely 
obviating the ability to recover correct yields. 

Readout laser power density. Figure 4C shows that 
variation of the laser power density within the range 0.2- 
20 kW/cm^ qualitatively has the same effect as changing the 
frametime. This can be easily understood, as, at low laser power, 
reduced SNR results in yield's overestimation, whereas at high laser 
power fluorescence traces get shorter due to increased photoex- 
citation, giving rise to insufficient time resolution. 

Molecular density. In contrast to in vitro samples, where a 
sufficiently sparse distribution of immobilized fluorescent mole- 
cules can be achieved to prevent overlap of fluorescence traces, in 
biological samples fluorescent labels typically decorate a fine 
structure of interest and their local density can be very high. This 
translates into a substantial density of active molecules per frame. 
Figure 4D shows that extraction of phototransformation yields 



remains successful at densities up to 7 molecules/frame/nm^ on 
our microtubule mimicking pattern (note that such densities on an 
essentially 1-D pattern cannot be directly compared to densities 
extracted from 2-D patterns [57]). At very high densities, an 
overestimation of the yields is observed, which results from 
rejection of strongly overlapped spots and heavy mix up of a large 
number of fluorescence traces (Figure 4D, inset). However, the 
results suggest that typical expression levels in biological cells do 
not need to be reduced to retrieve correct photophysical 
parameters, at least under weak 405 nm laser illumination. 

Number of frames. Figure 4D shows that retrieval of the 
correct photoconversion brightness necessitates a substantial 
fraction of the total number of molecules to be photoconverted 
during the experiment (>30%). This is explained by inspecting the 
exponential dependence of jV^i(Z) = jVo[l — exp(— ot^a/)] (see also 
Figure SIO). As long as NJJ) can be approximated by a linear 
regime MA{t)~Mo'X<pJ {(X.(pJ«i, (X is a constant) the unknown 
parameters jVo and p„ cannot be uncoupled and unreliable values 
are extracted. Finally, Figure 4F simply shows that for a fixed 
sample and in the presence of a heterogeneous laser beam a 
substantial number of activated molecules need to be recorded so 
that sufficient statistics allow reliable extraction of the photo- 
transformation yields. 

Overall, the analysis of simulated datasets presented above 
demonstrates that our procedure estimates the phototransforma- 
tion yields with suitable accuracy and robustness in a large range 
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Figure 6. Phototransformation yields of Dendra2-actin ex- 
pressed in mammalian cells and purified Dendra2 embedded 
into PVA film. Error bars indicate standard deviations. Number of 
samples per experiment n = 3-5. Cell samples are from different batches 
prepared independently. 
doi:10.1371/journal.pone.0098362.g006 

of experimental parameters. This constituted the starting point for 
choosing appropriate settings for real experiments and estimating 
the suitability of experimental data sets for analysis. 

Experimental data 

In vitro data. The phototransformation yields of Dendra2 
were first measured in vitro in PVA, under conditions of typical 
PALM data acquisition, with or without 405 nm activation (see 
Data acquisition section and Figure S2). The processing pipeline 
with the "tumbling" model was used. Results are reported in 
Figures 5 to 7. The obtained values are close to those that can be 
deduced from the phototransformation rates reported by Lee et al 
[8], although a different immobilization technique was used. With 
bleaching and blinldng-off quantum yields of (2.3±0.4) x 10~'^ and 
(8.6± 1.8)x 10~'\ respectively, Dendra2 exhibits a fair photosta- 
bility and has a relatively low tendency to blink, with on average 
one out of 3.7 molecules experiencing blinking. The rate of 
recovery from the dark state was measured to be 14±2 s 
corresponding to an average lifetime in the off state of 70 ms. Dark 
state recovery is essentially thermally driven and was found to be 
accelerated by 405 nm light only at laser intensities greater than 
(~10 W/cm^), which exceeds the values typically used in PALM 
acquisition, as previously demonstrated [8] (not shown). Likewise, 
blinking and bleaching quantum yields were not modified upon 



weak 405 nm laser light illumination. Upon treatment with 1,4- 
diazabicyclo[2.2.2]octane (DABCO), an antifading agent com- 
monly used in fluorescence microscopy to reduce photobleaching 
[64,65,66], the photobleaching and blinking-off yields slighdy 
decreased (04^,^: (1.5±0.5)xl0-^ <Pu,nk-off- (6.0± 1.4) xlO"''), 
while the blinking-on rate slightly increased [^Umk-m- 
18± 1.6 s~'), but all these changes were not significant (Figure 6). 

Finally, from the analysis of the cumulated number of 
photoconverted molecules over the time of PALM data acquisition 
(Figure 7), the photoconversion brightness of Deiidra2 was 
measured to be (0.9±0.2) X 10 ' M cm under 561 nm readout 
laser illumination. As expected, the brightness measured under 
405 nm illumination was substantially higher ((4. I±2.3)xl0~' 



) and, assuming an extinction coefficient of 10250 
at pH 7.4, as measured by ensemble spectroscopy, a 



M" 

M~-cm 

surprisingly low photoconversion quantum yield of 
(4.0±2.3)xl0 ^ could be deduced. Results are summarized in 
Table 1. Based on the coarse assumption of a wavelength 
independent quantum yield, the residual extinction coefficient of 
the neutral green state of the Dendra2 chromophore could be 
proposed to be ~ 1.0 M 'cm ' at 561 nm at pH 7.4. However, 
this statement would not hold if photoconversion in the anionic 
green state of Deiidra2 also occurs as receiitiy suggested [67]. 

Eukaryotic cells 

Next, we measured the phototransformation yields of Dendra2 
fused to actin in fixed mammalian Vero cells. A representative 
PALM rendered image is shown in Figure S12 and the results of 
analysis are shown on Figures 5 to 7 . The processing pipeline with 
the "fiKed-circular" model was used. Importaiidy, experimental 
conditions (notably signal-to-iioise ratio) were found to fall within 
the range suitable for accurate extraction of phototransformation 
yields (Figure 4). A photobleaching quantum yield of 
(2.6±0.4) X lO"'' was measured, close to the value obtained in 
PVA. In marked contrast, the blinking-off quantum yield 
(2.3±0.2)xlO ^ was found to largely exceed that measured in 
vitro, whereas the blinking-on rate (11.7±0.5 s~ ) was slightiy 
lower. Consequentiy, about half of the chromophores experienced 
blinking, a fraction significantly higher than in PVA, and also 
higher than the 30% ratio reported in Xenopus oocytes [21]. 

Strikingly, and in marked contrast with PVA, DABCO induced 
a significant enhancement of Dendra2 photostabUity in fixed Vero 
and HeLa cells. The photobleaching quantum yield was reduced 
nearly 2-fold ((1.4±0.2)xl0~'^), the blinking-ofiF quantum yield 
was decreased more than 3-fold ((6.6± 1.0) x 10 % whereas the 
blinking-on rate was increased by nearly 2-fold (22. 8± 1.7 s '). In 
fact, adding DABCO to the PBS medium used for cell imaging 
brought the phototransformation yields of Dendra2 close to those 
measured in PVA alone and nearly identical to those measured in 
PVA plus DABCO. However, DABCO did not cause an increase 
of mean photon output per molecule, possibly in agreement with 
the previously reported quenching effect of this substance [65] . 

Finally, the green-to-red in cellulo photoconversion brightness of 
Dendra2 was measured to be (0.9±0.4)xl0 ^ M cm and 
(1.4±0.6)xl0"' M"' cm , under 561 nm and 405 nm illumina- 
tion respectively (Figure 7). In the latter case, assuming an 
extinction coefficient of the neutral form of green Dendra2 at 
405 nm of 10250 M cm [36], a photoconversion yield of 
~(1.4±0.6) X 10~' could be deduced, even lower than in PVA 
(Table 1). 
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Figure 7. Representative cumulative photoconversion plots of purified Dendra2 embedded into PVA film (A, B) and Dendra2-actin 
expressed in mammalian cells (C, D); (illumination with 561 nm laser only (A,C) or in combination with 405 nm laser (B, D). Fitted 
curves are shown overlaid on experimental data. 
doi:1 0.1 371 /journal.pone.0098362.g007 



Discussion and Conclusions 

In this work, we have developed a reliable method to evaluate in 
cellulo the photo-transformation yields of fluorescent proteins used 
as markers in high-resolution localization microscopy. Through 
extensive simulations, we have determined the range of experi- 
mental parameters suitable for accurate determination of the 
photobleaching yield, the blinking-o£F yield, the blinking-on rate as 
weU as the photoconversion brightness of PTFPs. Tlie method was 
then applied to study the in cellulo photophysical behavior of 
Dendra2 and the results were compared to those obtained in PVA. 
An important conclusion is that the photophysical behaviors of 
fluorescent proteins measured in vitro and in cellulo may differ. 

Table 1. Photoconversion parameters of Dendra2. 



Quantitative localization microscopy therefore would benefit from 
a preliminary investigation of the fluorescent marker's bleaching 
and blinking behavior in biologically relevant experimental 
conditions. In the case of Dendra2 fused to P-actin in fixed 
mammalian cells at physiological pH, we found a bleaching yield 
of (2.6±0.4)xl0"\ a blinking-off yield of (2.3±0.2) xlO"^ and a 
thermally activated blinking-on rate of 11.7±0.5 s~'. Interesting- 
ly, when the antifading agent DABCO was added at 10 mM to the 
PBS-based medium in which mammalian cells were mounted, the 
photostabUity of Dendra2 was largely enhanced. Reversible 
blinking was decreased more than 3-fold and nonreversible 
photobleaching was reduced by about 2-fold. Furthermore, the 
average duration of Dendra2 transient dark states was reduced 







PVA film 


Mammalian cells 


Photoconversion brightness (561 nm) [M ^cm ^] 


(0.9±0.2)xl0"^ 


(0.9±0.4)x10"= 


Photoconversion brightness (405 nm) [M 'cm 


(4.1±2.3)xl0"' 


(1.4±0.6)xl0"' 


Photoconversion yield (405 nm)^* 


(4.0±2.3)xl0"^ 


(1.4±0.6)xl0"^ 



Mean values ± SD are shown. Number of samples per experiment n = 3. 

^: The high SDs are attributed to inhomogeneity and weakness of the 405 nm laser beam. 

*: Assuming an extinction coefficient <':405= 10250 MoP^ cm^V 

doi:10.1371/journal.pone.0098362.t001 
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from ~85 ms to ~45 ms. Thus, because less molecules wUl blink, 
and those that still blink will more rapidly return from their dark 
state, the use of DABCO in combination with Dendra2 may 
facilitate accurate molecular counting in quantitative PALM. 

While redox agents and oxygen scavenging systems have been 
carefully evaluated in STORM/ dSTORM localization microsco- 
py to ensure optimal blinking of the organic labels [47], the 
chromophores of fluorescent proteins are often considered 
insensitive to the local emironmcnt and FP-labeled samples are 
commonly prepared without any additives. However, recent data 
clearly demonstrated that the photophysical properties of fluores- 
cent proteins are sensitive to the environment [50,52,68,69]. Our 
results go along this line. 

Two main scenarios can be invoked to account for the 
observations reported in Figure 6. Firstly, DABCO is known to 
retard photobleaching of various fluorophores in immunofluores- 
cence samples, which was explained by hampering diffusion of 
oxygen and other reactive oxygen species (ROS) capable of 
destroying the fluorophore [64,65,66]. In addition, PVA is known 
to also hamper diffusion of oxygen [4] . Thus, a reduced oxygen 
level could possibly explain the enhanced photostabilit)' in PVA or 
in DABCO-enriched PBS if blinking and bleaching mechanisms 
under PALM experimentEil conditions are oxygen dependent. 
Alternatively, DABCO is also an electron acceptor and it has been 
shown that PVA may be involved in electron-transfer processes 
to/from organic dyes in their excited states, thereby affecting the 
blinking properties of these dyes [70,7 1] . Thus, our results could as 
well be explained by redox effects. 

Our previous investigations of IrisFP, a PTFP which shares high 
sequence-homology and structural-similarity with Dendra2 
[42,43,44], lead us to suggest that this latter scenario predomi- 
nates. Under high-intensity illumination as t>'pii:ally used in 
localization microscopy (>0.1 kW/cm^) structural studies indeed 
revealed that IrisFP photobleaching essentially follows a redox- 
driven process associated with the decarboxylation of the 
conserved Glu212 residue (IrisFP numbering) [44]. This photo- 
destruction pathway leads to a non-reversible 2-electron photore- 
duction of the chromophore, likely to be modulated by the local 
redox environment. The process is oxygen-independent, in 
contrast to the photobleaching pathway observed at lower 
illumination levels typical of standard widefield microscopy. 
Likewise, our investigations of blinking in IrisFP [42,43] silso 
suggested a redox process in which the FP chromophore 
undergoes a transient loss of electron conjugation at the methylene 
bridge, followed by proton transfer. DABCO and PVA could play 
a role in such mechanisms. Further investigations will need to be 
carried out to confirm or infirm this hypothesis. Preliminary in 
celluh experiments performed under anaerobic conditions suggest 
that indeed the measured Dendra2 phototransformation yields do 
not depend on oxygen significantiy. 

More generally, it will be interesting in the future to carry out 
more systematic studies of the effects of DABCO and other 
antifading agents such as e.g. trolox, on photoblinking and 
photobleaching of various fluorescent proteins. In addition, in the 
present study we have not looked at a possible dependence on 
DABCO of the Dendra2 photoconversion brightness. It wiU be 
informative in future experiments to investigate the possible redox 
dependence of green-to-red photoconversion in various PCFPs, as 
reported for example in mEos2 [51]. 

The phototransformation yields in Dendra2 reported in this 
work fall within the range of values determined for other PTFPs by 
various ensemble or single molecule methods. Although yields 
measured in the exact conditions of biological PALM experiments 
are expected to be most informative, their accuracy could be 



affected by various complications not accounted for in the present 
work, such as local heterogeneity of the laser beam, nonuniform 
background, highly anisotropic distribution of PTFP's dipole 
orientations, or misestimation of laser power density due to optical 
aberrations. Importantly, extinction coefficients are not measured 
in vivo, are pH dependent and are difficult to estimate accurately. 
Reliable measurements of phototransformation quantum yields 
require finding a good compromise between high signal-to-noise 
ratio of the recorded fluorescence and sufficient time-resolution to 
faithfully reconstruct single-molecule traces. In particular, we have 
shown that the limited time-resolution of PALM experiments may 
significantiy bias the retrieval of bhnking-off quantum yields, 
although such bias could possibly be corrected by generating 
calibration curves from simulated data. For example, in the case of 
the fixed Vero cells, applying the calibration curve of Figure 3 
(inset) would suggest a blinking-oflf quantum yields of 
(6.7±0.6)xl0~''' and (1.5±0.1)xl0~'' in the absence and presence 
of DABCO, respectively, meaning that DABCO induces a drop in 
blinking by a factor ~5. However, one should bear in mind that a 
proper calibration needs to faithfully account for a high number of 
experimental parameters, and we propose- that more thorough 
investigations should be conducted before- corrections can be 
confidently applied. Perhaps even more important, i:ell to cell 
variability always remains an issue to take into account and it 
cannot be excluded that sample preparation protocols as well as 
specific experimental conditions (e.g. temperature or presence of 
vitamins in the cell culture medium [52]) may exert an influence 
on the extracted phototransformation yields. Therefore well- 
defined artificial environments with single-molecule or even bulk 
measurements will continue to provide valuable reference values 
that should not be overlooked. 

Within all these limitations, our extensive simulations suggest 
that accurate phototransformation quantum yields can be 
retrieved if they fall within reasonable ranges and that reliable 
comparative studies between different samples or various envi- 
ronmental conditions can be conducted. 

Within the precision of our measurements, Dendra2 molecules 
displayed a homogeneous photophysical behavior, as NiUadi, ^on 
and }^'„jf histograms could be satisfactorily fit with mono- 
exponential models. This is in contrast with other findings made 
in vitro where recoverv' to the fluorescent state exhibited a biphasic 
behavior [8], the mechanistic interpretation of which remains 
unclear. Our data suggest that various conformational states of 
Dendra2, if they exist, interconvert at rates that exceed the 
measured macroscopic rates. 

Our simplified photophysical scheme (Figure 1) may also be 
refuu-d in the future to extract more accurate phototransformation 
yields. With more statistically significant data, multiphasic 
bleaching and blinking processes could be pinpointed [72], and 
the possibilit)' that bleaching may directly follow from blinking 
[59,60] could be investigated. One observation is that the 
reccn'ered green-to-red photocom'ersion brightness (or yields) for 
Dendra2 appears lower than expected based on the few reported 
bulk experiments on this protein [67] and other PCFPs [18,19]. 
The photoconversion quantum yields measured here under 
405 nm activation are of the same order of magnitude as the 
photobleaching or blinking-off quantum yields. This caUs for an 
experimental strategy in which bleaching and blinking in the green 
state of the fluorophores are taken into account, as weU as further 
photophysical processes that may be responsible for the generaUy 
observed limited overaU photoconversion efficiency, as recentiy 
shown [21]. 

Finally, the present study opens the door to a more systematic 
investigation of the influence of environmental parameters on the 
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photophysical behavior of markers or interest for single-molecule 
localization microscopy. The technique may for example be 
extended to the case of organic dyes [73], allowing a precise 
assessment of the effects of various blinking cocktails used in 
techniques such as dSTORM. 

Supporting Information 

Figure SI View of our Matlab-based PALM simulation 
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